Density and current response functions in strongly disordered electron systems: 
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We study consequences of gauge invariance and charge conservation of an electron gas in a strong 
random potential perturbed by a weak electromagnetic field. We use quantum equations of motion 
and Ward identities for one- and two-particle averaged Green functions to establish exact relations 
between density and current response functions. In particular we find precise conditions under which 
we can extract the current-current correlation function from the density-density correlation function 
and vice versa. We use these results in two different ways to extend validity of a formula relating the 
density response function and the electrical conductivity from semiclassical equilibrium to quantum 
nonequilibrium systems. Finally we introduce quantum diffusion via a response function relating 
the averaged current with the negative gradient of the averaged charge density. With the aid of 
this response function we derive a quantum version of the Einstein relation and prove the existence 
of the diffusion pole in the zero-temperature electron-hole correlation function with the long-range 
spatial fluctuations controlled by the static diffusion constant. 
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I. INTRODUCTION 

Low-energy physics of equilibrium systems with 
weakly interacting electrons is well understood both 
qualitatively and quantitatively. The relevant informa- 
tion about the equilibrium system is contained in one- 
and two-particle Green functions, in particular in their 
behavior around the Fermi energy. A number of re- 
liable approximate methods have been developed for 
the calculation of these functions. Among them, sys- 
tematic renormalized perturbation expansions based on 
the many-body Feynman diagrammatic technique have 
proved most effective. When impurities or quenched con- 
figurational randomness are added we still are able to 
describe equilibrium properties of such a system quite 
reliably within diagrammatic schemes and a mean-field- 
type coherent-potential approximation (CPA)ii 

The situation gets less straightforward if we come out 
of equilibrium. This is the case when the system is dis- 
turbed by a time-dependent external force but does not 
manage to reach a new equilibrium within the relaxation 
times of the experimental setup. We then have to de- 
termine the response of the system to the external per- 
turbation to obtain experimentally relevant data. How- 
ever, we do not have at our disposal many established 
methods to calculate nonequilibrium response functions. 
When the perturbation is weak, which is the most com- 
mon situation in practice, we can use linear response 
theory to calculate effectively the response functions.^ A 
fundamental tool for the calculation of response functions 
within the linear-response theory is the Kubo formalism!^ 
It determines how the response functions of the disturbed 
nonequilibrium system can be expressed in terms of two- 
particle Green (correlation) functions that are character- 
istics of an equilibrium state. Kubo formulas are means 
for the description of nonequilibrium systems with equi- 



librium functions. Unlike the response functions, the 
equilibrium Green functions obey equations of motion 
and are suitable for developing systematic approxima- 
tions within many-body diagrammatic techniques. 

Kubo formalism, however, provides different represen- 
tations for different response functions that we have to 
approximate separately. The problem that emerges with 
independent approximations for different response func- 
tions via different Kubo formulas is that we are usually 
unable to keep exact relations between response func- 
tions that may hold due to special symmetries of the 
system under investigation. Most pronounced case with 
"hidden" symmetries is the response to a weak electro- 
magnetic perturbation where we have to satisfy gauge 
invariance and charge and current conservations. Instan- 
taneous field-dependent deviations of charge and current 
densities from their equilibrium values are described by 
density and current response functions. Within linear- 
response theory these functions are obtained from inde- 
pendent Kubo formulas with density-density and current- 
current correlation functions, respectively. Although 
both correlation functions are derived from the same two- 
particle Green function, they are normally approximated 
independently according to the purpose they serve to. 

When quantum coherence is negligible we can calcu- 
late the transport properties from the coherent-potential 
approximation^ There is no contribution to the homo- 
geneous current-current correlation function beyond the 
single electron-hole bubble in the single-band coherent- 
potential approximation.** Hence most theories beyond 
the CPA, either on model or realistic level, use the Kubo 
formula for the electrical conductivity with the current- 
current correlation function to determine transport prop- 
erties of disordered solids 

On the other hand, when quantum coherence ef- 
fects are substantial and we expect the Anderson metal- 
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insulator transition, electrical conductivity is usually cal- 
culated from the electron-hole correlation function with 
the aid of the diffusion constant controlling its low-energy 
behavior near the diffusion poleA^ The diffusion pole 
is crucial and enables application of scaling arguments 
and the renormalization group approach to Anderson 
localization.^ 

There is a number of more or less heuristic arguments 
in the literature that relate the density response with the 
conductivity.^*^ They are essentially based on gauge in- 
variance and charge conservation of a semiclassical equi- 
librium description of an electron gas exposed to an elec- 
tromagnetic field. In weakly disordered quantum systems 
(described by continuum models) a relation between the 
density response and conductivity can formally be de- 
rived as followsiiS Gauge invariance is used to relate the 
external scalar potential with the electric field E = — 
The current density generated by the external field then 
is 

j(q, uj) = o-(q, uj) ■ E(q, uj) = -icr(q, w) • q (^(q, uj) (1) 

where cr(q, w) denotes tensor of the electrical conduc- 
tivity. Charge conservation is expressed by a continuity 
equation. In equilibrium we can use its operator form 
that follows from the Heisenberg equations of motion for 
the current and density operators. For Hamiltonians with 
quadratic dispersion relation we have 

e9tn(x,i)+V-J(x,i)=0. (2) 

Energy-momentum representation of the continuity 
equation in the ground-state solution is 

— ia;e(5n(q, Lj) + iq • j(q, cj) = . (3) 

We have to use a density variation of the equilibrium 
density, i. e., the externally induced density 6n{q,uj) — 
n(q, — no in the continuity equation with averaged 
values of operators. From the above equations and for 
linear response Sn{q,uj) = — ex(q, w)(y9(q, w) we obtain 
in the isotropic case 

-ie'^uj 

c^(q,^^) = — 5— x(q,^ ■ (4) 
q 

The derived equality formally holds for complex func- 
tions without restrictions on momenta or frequencies. 
Frequencies can, in principle, be even complex. Rela- 
tion Q is often taken as granted for the whole range 
of the disorder strength and used for the definition of 
the zero-temperature conductivity when describing An- 
derson localization transition. ^'^ 

Although the above derivation may seem very general 
it suffers from a few fiaws. First, the operator continu- 
ity equation Q cannot be directly used out of equilib- 
rium. A new term due to causality of the response func- 
tions enters the continuity equation when the equilib- 
rium is disturbed. Moreover the noncquilibrium density 
and current operators in linear-response theory no longer 



obey Heisenberg equations of motion with the perturbed 
Hamiltonian. The perturbation is decoupled from the 
equilibrium Hamiltonian and is treated only to linear or- 
der. Second, strongly disordered electron systems can 
be desribed only by lattice models with a nonquadratic 
dispersion relation. Continuity equation ^ is then to 
be modified even in equilibrium. Hence the above, in 
the literature broadly disseminated reasoning leading to 
Eq. cannot be fully trusted in quantitative studies of 
strongly disordered systems out of equilibrium. 

Presently there is no reliable theory of strongly disor- 
dered electrons beyond the mean-field CPA. Since mean- 
field approximations do not include vertex corrections to 
the electrical conductivity they are not suitable for the 
investigation of localization effects in three spatial di- 
mensions. Unlike low dimensions (d < 2), Anderson lo- 
calization in d = 3 may occur only in strongly disordered 
systems. Anderson localization in bulk systems has not 
yet been understood or quantitatively described in a sat- 
isfactory manner. Its quantitative description demands 
bridging the gap between the mean-field (CPA) trans- 
port theory and theories for weakly and strongly local- 
ized electron states. It is clear that such an interpolating 
scheme should be based on advanced approximations for 
two-particle functions. 

Recently a diagrammatic method for summations of 
classes of two-particle diagrams has been proposed. ^"'^ 
This theory has a potential to interpolate between the 
mean-field and localization theories provided exact rela- 
tions between the density and current response functions 
have been proved for this scheme in the whole range of 
the disorder strength. Criteria for validity of relations 
between density-density and current-current correlation 
functions in approximate treatments of strongly disor- 
dered systems have not yet been established. When we 
want to go beyond the CPA we still have to establish such 
relations in the metallic regime of tight-binding models 
with extended electron states. Simultaneously we have 
to formulate conditions under which the derived relations 
hold or may be broken. 

The aim of this paper is threefold. First, we de- 
rive various relations between current-current, current- 
density, and density-density correlation functions for 
strongly disordered electrons out of equilibrium, i. e., be- 
yond the reach of the existing derivations. Second, we 
articulate conditions under which these relations hold or 
to which extent they may be broken in quantitative treat- 
ments and show when the conductivity can be calculated 
from the density response function and vice versa. Third, 
we introduce a generalization of diffusion via a quantum 
response function, test validity of the Einstein relation 
for this quantum diffusion, and relate quantum diffusion 
with its classical counterpart and with the diffusion pole 
in the electron-hole correlation function. 

To reach this goal we use Kubo formalism and av- 
eraged many-body Green functions. In this approach a 
weaker form of the continuity equation expressed in terms 
of one- and two-particle Green functions replaces the op- 
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erator identity. It is derived from equations of motion 
for Green functions and Ward identities between them. 
We use two Ward identities that are not fully equivalent. 
One expresses conservation of probability and the other 
reflects charge conservation. Each Ward identity is used 
in a different manner to relate the conductivity with the 
density response. 

The layout of the paper is as follows. In Sec. ^ we 
summarize the definitions and useful representations of 
the density and current response functions. In Sec. IIIII 
we show how the current-current correlation function 
emerges from a momentum expansion of the density 
response function. Then in Sec. II VI we derive continuity 
equations expressed in Green functions and use them to 
show how the density response function can be revealed 
from the conductivity. Based on the exactly derived re- 
lation between the density and current response function 
we introduce in Sec. quantum diffusion, derive the 
diffusion pole in the hydrodynamic limit of the density 
response function and relate the diffusion constant to the 
conductivity. In Sec. IVII we illustrate general results on 
an exactly solvable limit of infinite spatial dimensions. 
In Appendix we discuss assumptions for and the range of 
validity of the VoUhardt-Wolfie-Ward identity used in the 
derivation of the continuity equation for Green functions. 



II. DENSITY AND CURRENT RESPONSE 
FUNCTIONS IN DISORDERED SYSTEMS 

A simplest description of the electron motion in im- 
pure, weakly correlated metals is provided by a tight- 
binding Anderson model. It assumes noninteracting spin- 
less electrons moving in a random, site-diagonal potential 
Vi. Its Hamiltonian reads 



Had = ^ t,j4cj -H^F.cIq. 



(5) 



The values of the random potential Vi are site- 
independent and obey a disorder distribution p{V). I. e.. 



a function depending on the random potential Vi is av- 
eraged as 



dVp{Vi}X{Vi} . 



(6) 



The averaged two-particle propagator (resolvent) is de- 
fined as the averaged product of one-particle propagators 



zil-t-V 



Z2\-t-V 



(7) 



Fourier transform to momentum space is not uniquely de- 
fined, since due to momentum conservation we have only 
three independent momenta. For our purposes we choose 
the following notation and definition of the Fourier trans- 
form 



ijkl 



The density response function is determined from 
the following thermodynamic representation of the Kubo 
formula— 



- ^ ksT ^ g'^I, {iujn, iujn + iv,n] q) • (9) 



kk' 



Here LOn = (2n -f 1)t:T are fermionic and Vm = 2rmrT 
bosonic Matsubara frequencies at temperature T. Due 
to the analytic properties of the two-particle Green func- 
tion we can perform analytic continuation from imagi- 
nary Matsubara to real frequencies using contour inte- 
grals. Only the contribution along the real axis remains 
and we obtain 



X(q, -f zO+) ^—Y^j^ — {[f{E + f{E)] G^^iE, + c.; q) 



f{E)G^^,{E, E + oj-c^)- f{E + Lo)G^C'{E, E + uj;c^)] 
I 



(10) 



We used the following definition 

Gt5{E, E + u-ci)^ G^^l {E - *0+, E 



zO+;q) 



and analogously for the functions G^^ and G"^"^, where 
both energies are from the upper and lower complex 



energy half-plane, respectively. We denoted the Fermi 
function/(£') = [1 + exp{/3(£^ — /i)}]^^ with the chemical 
potential p,. 

The two-particle Green function is generally deter- 
mined from a Bethe-Salpeter equation in which the input 
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is a two-particle irreducible vertex. The latter is known 
only approximately, except for special limits. However, a 
few specific elements of the full two-particle Green func- 
tion can be evaluated without knowing the two-particle 
irreducible vertex. It is sufficient to know only the one- 
particle Green function to determine them. Determina- 
tion of two-particle functions from one-particle ones is 
enabled by Ward identities. The Ward identity relating 
the averaged one- and two-particle Green functions reads 

^ E G^^l (^1' 0) = [G(k, zi) 

-G(k,Z2)]. (11) 

and was proved for the first time within the coherent- 
potential approximation by Velicky^ It is a nonpertur- 
bative identity valid quite generally beyond the CPAiii 
As a consequence of this identity we obtain vanishing of 
the density response function for a homogeneous global 
perturbation, i. e., for q = 0. Using lfTT|l in ifTUIl we easily 
find 

X(0,w + z0+) 

- G-'iB + u)] - /(E) [G-'iE) - G'iE)] } = . (12) 

Another situation where we do not need to know the 
two-particle irreducible vertex is the static limit, oj = 
0. First, the static density response function is real and 
reads 

^(q' ^)-^T. —fiE)^GS^;{E, E- q) (13a) 
kk' 

that in the limit q ^ goes over to 



with np as the density of states at the Fermi level and 
G^(cj) = iV-iX;qG^(q-'^)- Note that Eqs. ^ hold 
only in the limit Lu/q ^ 0. In the inverse case, q/to — > 0, 
Eq. (|12|l applies. Noncommutativity of the limits u 
and (J — > indicates that the point q = 0, uj — is 
not analytic in the same manner as in the Fermi liquid 
theoryii^ 

The current-current correlation function 11^^ (q, t) re- 
lates the current of the system j(q, i) perturbed by an 
external vector potential A(q, t). Tensor of the elec- 
trical conductivity (7^/3 (Qi t) determines the current re- 
sponse to an external electric field. It can be obtained 
from the current response function by using a relation 
E(x, f) = A(x, t). In the Fourier representation we 
obtaini^ 



cr„/3(q,t^+) = - [nc/3(q,t^+) - nQ^(q,0)] (14) 

UJ 

with = uj + iO~^ . The compensation term 0) on 

the right-hand side of Eq. 114|) is real and was introduced 
to warrant finiteness of the complex conductivity. This 
additive term is not generated from the Kubo formula 
with a commutator of current operators but originates 
from gauge invariance of the electron gas in an external 
electromagnetic field»i& According to Eq. (|14|l there is no 
current in the system with a static vector potential and 
the static (complex) conductivity can be defined only via 
a dynamical one from the limit u —> 0. 

Unlike the density response function the current re- 
sponse function cannot be simplified in any limit and 
to determine the conductivity from Eq. (|14l) we have to 
know the two-particle irreducible vertex. The current re- 
sponse function alike the density response function is a 
complex quantity with the same analytic properties. Of 
interest for us are the parallel components that for real 
frequencies read 



a«„(q,c^) = -__ ^ [z;„(k + q/2) - v„{-k + q/2)] [w„(k' + q/2) - «„(-k' + q/2)] 

k,k' 

^ ^ { [/(i^ + ^) - f{E)]G^3 {E. + c.; q) + f{E)G^S' (E, E + co; q)) 

~f{E + u;)G^^,{E,E + oj;q)-f{E)[G^^iE,E;<i)-G^^,{E,E;q)]} . (15) 
I 

We denoted the group velocity ^^(k) = m-^d/dkaeik), III. DYNAMICAL CONDUCTIVITY 

where e(k) is the dispersion relation of the underlying CALCULATED FROM THE DENSITY 

lattice. RESPONSE 



As first we show how the dynamical conductivity can 
be obtained from the long-range behavior of the density 
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response function (hydrodynamic limit of small transfer 
momentum). To this purpose we use the Velicky-Ward 
identity, Eq. Hll|) . It enables us to evaluate averaged ma- 
trix elements of the two-particle Green function with zero 
transfer momentum in terms of the one-particle propa- 
gator. Zero transfer momentum is a severe restriction 
on the applicability and utilization of this Ward identity. 
Physically we are interested in the hydrodynamic limit, 
i. e., the asymptotics q ^ 0. We are unable to extend 
identity Ijllfl beyond the homogeneous case. However, 
if we assume analyticity of the asymptotics q —> 0, we 
can use momentum q as an expansion parameter and 
investigate the hydrodynamic limit perturbatively. The 
hydrodynamic limit is analytic if frequency lu ^ 0, i. e., 
we are in the regime q/uj <^ 1. We then can expand the 
density response function in powers of momentum q. 

We define two configuration-dependent resolvent op- 
erators 



G±(z) 



(16) 



where Ait ± is a difference in the dispersion rela- 
tion of the two resolvents and is defined via its matrix 
elements (k|Aitlk') = S{k' - k) v(k) ■ q/2, (k|A2t|k') ^ 



S(k.' - k) (q • VA;)^e(k)/8. We wih be interested in the 
following function of two energies 



/(z+,z_;q) = -Tr[G+(z+)G_(z_) 



(17) 



that can be, after averaging over the configurations of the 
random potential, expressed via the two-particle Green 
function 

(/(z+,z_;q)),, = i^^G^'J,(z+,z_;q) . (18) 



kk' 



We assume that the expansion in the operators Ai^2t 
commutes with the configurational averaging. We first 
expand the quantity / and then average the series term 
by term. It is sufficient for our purposes to expand only to 
second order in momentum q. This precision determines 
the leading small-momentum behavior. 

Expanding quantity / in Ai,2t we have to keep the 
order of operators in the products, since Ai^2^ and the 
resolvent G{z) do not commute. The expansion to second 
order reads 



/2(z+, z_q) = Tr {g(z+)G(z_) + G(z+) [Ait + A2t\ G(z+)G(z^) 

- G(z+)G(z_) [Alt- Aat] G(z_) - G(z+)Ai?G(z+)G(z_)Ai?G(z_) 

+ G{z+) AitG{z+)AitG{z+)+G{z^)AitG{z^)A^t G(z_)} . (19) 

Each direct product of the resolvent operators G{z±)G{zz^) can be simplified using identity (|ll|l . Doing this conse- 
quently we end up with a sum of products of two resolvents. We have three different terms to analyze: 



Tr {g(z+) [{Ait + A2t) G(z_) - G(z+) (Aif- A2?)] G(2_)} 

{Tr [(g'(z+) + G'(z_)) Alt 



Tr 



G'(z+)-G'(z_)j A2t 
2 



z+ — z_ 



-Tr 



Ai 



f(G(z+)-G(z_))] I , (20a) 



TTlG(z+)AitG(z+)G(z)AitG(z\] = ^ (Tr \G(z+)AitG(z+)Ait 

I i (z+ — z_ )^ L L 



(Z+ -Z_)2 

+ Tr 



G{z-)ArtG{z-)Ait -2Tr G(z+)AitG(z_)Aif | , (20b) 



Tr{G(z+) Ai?G(z+)AifG(z+) + G(z_)AifG(z_)Aif G(z_)} 



|Tr G'(z+)AirG(z+)Aif - Tr G'(z_)AirG(z_) Ait 



z^ — z_ 



-Tr 



G(z+)-G(z_) Aif G(z+)-G(z-) Alt \ . (20c) 



We insert Eqs. H20|) in expansion 119|) and average over the configurations of the random potential Vi. To 
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recover the density response function we have to Hmit 
the complex energies to the real axis from above or be- 
low with which we distinguish the causality. We define 
correlation functions that can be represented as a trace 
of the averaged two-particle Green function 



-y 

kk' 



Gti,iE,E + u;;<D . (21) 



where a, b stand for A, R depending on whether the imag- 
inary part of the corresponding frequency argument is 
positive or negative. We then can represent the leading 
asymptotics in small momenta q of the averaged electron- 
hole correlation function 



<i>^«(q, to) - $^«(0, u^) = -^ [{Oe {G'^iE + u;) + G^'iE)) q • v)] 



1- [{dE {G^E + c) - G^{E)) (q . V)^6)] + \ [((G«(i? + c.) - G^{E)) q • v)] 
^ 5] q • v(k) q . v(k') {d, [G^^,{E + e, E; 0) - GS^,{E + lu, E + uj + e;0)] 



kk' 



+ I [G^^^{E,E;0) + G5^;{E + u;,E + u;-0)-2G^^,{E,E + u;;0)]'^ (22a) 



and analogously of the electron-electron [^"^^) and hole- hole ($ ) ones that can be written generically as 

$|»(q, u;) - $|»(0, ^) = We {G'^{E + c.) + G'^^E)) q • v)] 

+ ^ {{dE [G^iE + L0)- G'^iE)) (q . Vfe)] + ^ [{{G^iE + co) - G^iE)) q • v)] 
^ ^ q • v(k) q . v(k') {4 [G^k' (E + e, E; 0) - G^k' {E + to, E + uj + e;0)l^. 



kk' 



+ - [G^k' (E. E- 0) + Gil, {E + u,E + uj-0)- 2G£^, {E, E + uj-Q)] \ ■ (22b) 

UJ 

I 



In the above equations we used angular brackets to de- 
note summation over fermionic momenta in the first Bril- 
louin zone of one-electron functions, that is 



(23) 



Note that only the term with lo~'^ from the electron-hole 
correlation function diverges in the limit w — > while 
in the electron-electron (hole-hole) functions it remains 
finite even in the static limit. This is a manifestation 
of the diffusion pole missing the latter two correlation 
functions. 

In Eq. different powers of momentum q and fre- 
quency UJ appear. There is no restriction on validity of 
Eas. (|22l) in frequency but it holds only for small mo- 
menta, more precisely only perturbatively up to second 
order. Equations H22(l establish relations between the 
density-density and the current-current correlation func- 
tions in the asymptotic limit q = Q. Generally, however, 
these two functions are not directly proportional, since 
the right-hand sides of Eos. 1)22(1 contain one-particle con- 
tributions. They cancel each other if we combine the 



electron-hole with electron-electron and hole-hole corre- 
lation functions appropriately to build up the density re- 
sponse function. 

For isotropic situations we define in the hydrodynamic 
limit 



g{uj) = -ilim ^x(q,w) . 



(24a) 



This quantity can be generalized to anisotropic cases as 



-lUJ 



(24b) 



Function g{u!) measures leading long-range correlations 
of the density response function. Its real part can be iden- 
tified with the diffusive conductivity or mobility of the 
system. It is now easy to find an explicit representation 
of this function using Eqs. (|22|l . All contributions except 
for the last terms on the right-hand sides of Eqs. H22|l can- 
cel and we obtain equality aaa{^,i^) = e.'^daai^) being 
just Eq. Q in the limit q = 0. Note that in the course 
of the derivation we had to sum contributions from the 
electron-hole, electron-electron, and hole-hole correlation 
functions and integrate over energies. 
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If we resort to the low-frequency limit w — > 0, the dif- 
fusion pole in the electron-hole correlation function dom- 
inates in the density response and we recover the conduc- 
tivity directly from <i)'^^. It is interesting to note that in 
this static limit we derive the mobility (dc-conductivity) 
with contributions from G^^ and G^^ solely from the 
electron- hole correlation function 



IV. DENSITY RESPONSE FUNCTION 
CALCULATED FROM THE CONDUCTIVITY 

We showed in the preceding section how the conduc- 
tivity can be revealed from the small-momentum behav- 
ior of the density response function. In this section we 
express the current-current correlation function in terms 
of the density-density correlation function. The aid we 
use for this task is a continuity equation and another 
Ward identity 

- s(k_,z_) = -^^Akk'(^;+,^-;q) 

k' 

X [G(k;,z+)-G(k'_,z_)] . (25) 

proved for retarded and advanced functions by VoUhardt 
and Wolfle in Ref. IS We denoted k± = k ± q/2. 



Note that identity for q = 0, z+ = Ep + iO+, z_ = 
Ep — iO~^ was already used in Ref. [l^- Identity (|25ll is 
related to Eq. , however, the two Ward identities are 
identical neither in the derivation nor in the applicability 
and validity domains. The latter holds for nonzero trans- 
fer momentum q, i. e., for an inhomogeneous perturba- 
tion while the former only for q = 0. On the other hand, 
we show in the appendix that Eq. (|25|l . unlike Eq. 
holds only perturbatively within an expansion in powers 
of the random potential. I. e., we can prove this iden- 
tity for nonzero transfer momentum q only by assuming 
that perturbation expansions for the self-energy E and 
simultaneously for the two-particle irreducible vertex A 
converge. It was shown earlier that in the homogeneous 
case, q — 0, the VoUhardt- Wolfle- Ward identity follows 
from the Velicky-Ward oneiii 

To derive a continuity equation for Green functions 
we start with an equation of motion for the two-particle 
Green function. It is a Bethe-Salpeter equation where 
the input is a two-particle irreducible vertex. We have 
three possibilities (topologically distinct scattering chan- 
nels) how to construct a Bethe-Salpeter equationiii For 
our purposes the electron-hole channel is the relevant or 
most suitable one. There the Bethe-Salpeter equation in 
momentum representation reads 



GjfJ, (z+, z_ ; q) = G(k+ , z+)G(k_ , ) 



(26) 



where Akk' (2^+ , ; q) is the two-particle irreducible ver- 
tex from the electron-hole channel. One-electron propa- 
gators G(k±,z±) = [z±-E(k±q/2,z±)-e(k±q/2)]-i 
are the averaged resolvents. 

The product of the one-electron propagators from the 
right-hand side of Eq. H26|) can be decomposed into 

G(k+,z+)G(k_,z_) = 

AqG(k; z+,z-) 

Az - Aqi;(k;z+,z_) - Aqe(k) ^ ' 

where we denoted 

Aqe(k) = e(k+)-e(k_) , (28a) 



AqG(k; z+, z_) G(k+, z+) - G(k_, z_) . (28b) 



and analogously the difference AqE(k; z^, zJ). We mul- 
tiply both sides of the Bethe-Salpeter equation (|26|l by 
the denominator from the right-hand side of Eq. H27|l and 
assume validity of Eq. I|25|) relating the two-particle ir- 
reducible vertex A with the one-particle irreducible self- 
energy S. We then obtain a "difference" equation of 
motion 



[Aqe(k) - Az] G^'J,(z+, z_; q) = AqG(k; z+, z_)J(k - k') + ^ J] Akk" (2+, 2-; q) 



AqG(k;z+,z_)G^l)k'(^+,^-;q) -Gl,"k'(^+>^-;q)AqG(k";z+,z_) . (29) 



.(2) 
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We now define correlation functions generalized to com- 
plex frequencies describing density-density and density- 
current correlations by summing over the fermionic mo- 
menta in Eq. (|25J| 

$(zi,Z2;q) = jp^G^^l,{zi,Z2\q) , (30a) 

kk' 

$,(zi,Z2;q) = ^^Aqe(k)4k'(^i.^2;q) , (30b) 

kk' 

|.,(zi,Z2;q) = ^^4'J,(zi,Z2;q)Aqe(k') . (30c) 
kk' 

The contributions from the two-particle irreducible ver- 
tex A in Eq. (|29|l cancel each other provided the two- 
particle irreducible vertex is symmetric, i. e., Akk' = 
Ak'k- If so, we end up with a continuity equation relat- 
ing the generalized density-density and density-current 
correlation functions 

$e(2+,^-;q) - Az$(z+,z_;q) 

k 

Another continuity equation can be derived by multiply- 
ing Eq. H29I) with the energy difference and summing over 
the fermionic momenta. We obtain an equation relating 
the generalized current-current and density-current cor- 
relation functions 

iz+ , ; q) - Az^, (z+ , z_ ; q) 

k 

Combining the above two continuity equations we obtain 
the resulting relation between the generalized current- 
current and density-density correlation functions replac- 
ing the operator continuity equation 

$ee(2+,z-;q) - (Az)^$(z+,z_;q) 

= ^E^q^(l^;^+'^-)[^^ + ^q<l^)] • (31c) 

k 

Actually, correlation function (f>ee is strictly speak- 
ing not the current-current correlation function, since it 



is not the momentum or velocity that is summed with 
the two-particle Green function. We used an energy 
difference Aqe(k) from Eq. (|28a|l . Only in the case of 
quadratic dispersion relation we have Aqe(k) = q • k/rn 
and the energy difference is proportional to the group 
velocity. Otherwise equals the current-current corre- 
lation function only in the small-momentum limit, g — > 0. 
We use 

Aqe(k) « q • v(k) (32) 
to convert the energy difference in the hydrodynamic 
limit to a multiple of the group velocity v(k) and de- 
note 

$c«/3(2+,2;_;q) 

= ]^E""(kK(k')GLk'(^+>^-;q) • (33) 

kk' 

With the aid of Eq. (|33|l we rewrite continuity equation 
(pTTcll to 

J2 g"/<J>aM^+, q) (^"^^G{z+, z_)) 

a/3 a 

= Az [Az$(z+, z_; q) + (AqG(z+, z_))] (34) 

where the angular brackets stand for the summation over 
the fermionic momenta from the first Brillouin zone as 
defined in Eq. (|23|) . In the isotropic case we obtain 

<^aa{z+,z^;q) = ^q- (vAqG(z+, z_)) 

+ ^[Az<I>(z+,z_;q) + (AqG(z+,z_))] . (35) 
9 

It is the most general relation between the current- 
current and density-density correlation functions for ar- 
bitrary complex frequencies and non-zero momenta. It 
holds for finite momenta as far as the dispersion law re- 
mains quadratic. For general dispersion relations the va- 
lidity of Eq. (|35|l remains in the asymptotic regime of 
small momenta. 

To return to measurable quantities we limit the com- 
plex frequencies to the real axis. To derive exact relations 
we first define a generalization of the conductivity with 
the correlation functions 



o-«(q,t^) 



dE 
2ttuj 



{ [f{E + f{E)] $f/(£;, E + u;;<i) + f {E)<^>^/' {E , S + c.; q) 

-fiE + Lo)^i^{E,E + Lu;q)-fiE)[<P^J'iE,E;ci)~<fi^iE,E;q)]} . (36) 
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We use continuity equation (|31c|) to relate (Tee and x a-nd find 

„2 pca 7771 

aee(q,t^) +*e2u;x(q,c.) = - i^f{E) [^^^''[E^E-ci) - <^i^{E,E-c^)] 

r ^/(i^)^E^i^(k)[3G«(k+,£;)-3G«(k_,i?)] = 0. (37) 



Vanishing of the right-hand side in Eq. p7ll can be expUc- 
itly manifested when continuity equation (|31c|) is appUed 
to the two-particle correlation functions '^f^^ and 
The two terms on the right-hand side of Eq. (|37|l add to 
zero. 

We can extract the electrical conductivity from 
Cee (q, ^) in the limit q ^ or for quadratic dispersion 
law by replacing $ee = La la'^aa leading to (Tec(q, cj) = 
L^ g^crQ,Q.(q, a;). In the isotropic case we then reveal 
Eq. from Eq. (|37() . It is important to state that this 
proof of validity of Eq. Q is strongly based on the Ward 
identity (|25|l for nonzero momenta. The latter can be 
proved only perturbatively, see Appendix, and hence it 
is not clear whether Eq.l^J holds beyond the perturba- 
tive regime near and beyond the Anderson localization 
transition. 



V. DIFFUSION, DIFFUSION POLE, AND 
EINSTEIN RELATION IN QUANTUM SYSTEMS 

In the preceding sections we proved validity of Eq. 
for arbitrary frequencies and small momenta from the re- 
gion where quadratic dispersion relation is accurate in 
strongly disordered lattice systems. It is a useful rela- 
tion in particular in approximate theories. It enables one 
to approximate only either the density response func- 
tion or the conductivity and to determine the other one 
from Eq. Q. However, both quantities in Eq. are 
response functions that are not well suited for system- 
atic approximations. We would prefer exact relations or 
identities for equilibrium correlation or Green functions. 
This can be achieved when we introduce diffusion. Dif- 
fusion in classical systems is defined by Pick's law relat- 
ing the current with the negative gradient of the charge 
densityiSii^ The diffusion constant introduced in this way, 
i. e. via Pick's law, has in classical physics an intu- 
itive phenomenological character without a proper micro- 
scopic backing that should come from quantum physics. 
The study of diffusion in quantum disordered systems 
was launched by the seminal work of Anderson analyzing 
destructive effects of quantum coherence on diffusion of 
disordered electronsiSS Since then diffusion has become 
part of transport studies of disordered and interacting 
quantum itinerant systems.---' ^•^'^-'^ However even in these 
quantum approaches diffusion is introduced either via 
(self-consistent) perturbation expansions leading to the 
diffusion pole in the electron-hole correlation function or 



via the semiclassical limit where the electron-hole corre- 
lation function becomes a Green function of the classical 
diffusion equationiiii 

Here we propose a nonperturbative way to define dif- 
fusion from first principles. We introduce a quantum gen- 
eralization of Pick's law with a diffusion response function 
relating the averaged current density with the negative 
gradient of the charge density of a perturbed system. We 
hence use 

(j(x,t)}^^ ^-e r dt' I d^x'D{^~^\t-t') 

J —00 J 

X V(5n(x',i'))a. (38a) 

as a definition of the diffusion response function D{x,t) 
being an extension of the diffusion constant to which it 
should reduce in the semiclassical limit. In the above def- 
inition we used a perturbed charge density (5n(x, t) that 
is not the total change of the charge density 5n(x, t) = 
n(x, t) — uq. It is a "dynamical" density change obeying 
a condition 

/oo 
dt{6n{x,t))av^0 (38b) 
-oc 

saying that the total change at any place during the entire 
time evolution is zero. This restriction is dictated by the 
same condition obeyed by the averaged current density 
defined from the conductivity (|14|l . 

Definition of diffusion (|38|l is nonperturbative, has 
quantum character and does not demand any assump- 
tion about the behavior in the semiclassical limit or the 
existence of the diffusion pole in the electron-hole corre- 
lation function. Moreover, it offers a microscopic deter- 
mination of the diffusion function from first principles, 
i. e., quantum response functions of the perturbed sys- 
tem and via the Kubo formalism from correlation and 
Green functions of an equilibrium solution. It is hence 
directly related to measurable quantities. 

In Pourier representation we obtain for the averaged 
current 

(j(q, ^)av = -«eq£'(q, uj) {Sn{q, uj))^^ . (39) 

The dynamical charge density Sn is induced by an exter- 
nal dynamical scalar potential, which can be expressed 
in Pourier representation as 

(<5n(q, u))^^ = - e [x(q, u;) - x(q, 0)] ^^(q, u;) . (40) 
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It is easy to verify that with this definition Eq. 138b|l is 
fulfilled. We use Eq. ((TJ together with Eqs. and (gDJ 
to find an identity between the conductivity, diffusion 
and density response 



(41) 



We have derived a general relation between quantum 
diffusion and the electrical conductivity. It is now to 
show that quantum diffusion response function D{q,Lu) 
behaves in the semiclassical limit as expected from the 
diffusion constant in the diffusion equation. Namely, 
we have to show that the homogeneous static diffusion 
D{0, 0) obeys Einstein's relation and enters the diffusion 
equation in the semiclassical limit. To come to it we 
use Eq. Q to exclude the density response x(q, i^) from 
Eq. (|41(l and Eq. H13a|l to determine x(q, 0). Doing so 
we obtain in the isotropic case 



cr(q,w) 

/°° rIF' 
— /(£;)SJ$i^(q,0) . (42) 
-oo 7^ 



If we define the homogeneous dynamical conductivity 
(T{io) = a(0,uj) (in the same way also the dynamical dif- 
fusion) and use the Velicky-Ward identity (|ll|l in the 
limiting case zi — Z2 — we obtain a dynamical gener- 
alization of the well known Einstein relation 



-oo 



that at zero temperature reduces to 
(T(aj) — e^npD{uj) 



(43b) 



A classical version of this formula was proved by Ein- 
stein for the Brownian motion of a particle in a random 
medium^S and later on it was re-derived in the frame- 
work of nonequilibrium statistical mechanics by Kubo.'^ 
Here we derived Einstein's relation for quantum response 
functions and showed that it is a consequence of gauge 
invariance of the system, namely Eq. 

It is not yet clear from the above reasoning how the 
diffusion function £'(q, lj) is related to the classical dif- 
fusion equation and to the diffusion pole. The existence 
of the diffusion pole in quantum systems is usually de- 
duced from the semiclassical limit leading to the diffusion 
equation. We can prove the existence of the diffusion pole 
entirely from first quantum principles and from the quan- 
tum generalization of Pick's law, Eq. (|38|) . We use Eq. Q 
to exclude now conductivity (7(q, lj) from Eq. 1)41(1 . It is 
easy to come to a representation of the density response 
function 



ZJ(q,c^)x(q,0)g^ 
^(•^'"^ = ~ — TTTi 



(44) 



known from other treatments of diffusionjiS Eq. 144(1 
holds for arbitrary frequencies and momenta within the 
range of quadratic dispersion relation. This represen- 
tation directly indicates a singularity in the density re- 
sponse function in the limit uj 0, q 0. Since the 
order of the limits is relevant for the result, we have to 
specify it explicitly. To single out the singular contribu- 
tion in the density response function we have to choose 
uj/q ^ 1, that is, the inverse situation to Sec. Illll and to 
Eqs. (PI). 

The singularity in the density response function leads 
to a pole in the two-particle Green function. It is demon- 
strated in the low-energy and small-momentum asymp- 
totics of the density response function at zero tempera- 
ture. We find from Eq. ifTUI) 



X(q, CO) = x(q, 0) - ($^^(q, 0) + 0{q")) 

+ 0{u;^). (45) 

Using this asymptotics in Eq. I(44() we obtain an ex- 
plicit manifestation of the diffusion pole in the zero- 
temperature electron-hole correlation function 



i>^?(q,^) 



27m i? 



'iuj + Dq"^ 



(46) 



representing the singular part of the electron- hole corre- 
lation function. We denoted a diffusion constant D = 
limg^o D{q, 0). This representation was derived for zero 
temperature in the asymptotic limit g —> 0, w — > with 
the restriction to/q <^ 1. To extend it to the opposite 
limit, q/uj <C 1, we have to show that the diffusion func- 
tion i?(q, w) is analytic in the limit g — > 0,a; — > 0. It is 
not a priori clear that the diffusion constant in Eq. I(4ti() 
equals the diffusion constant obtained from the inverse 
order of limits D = limi^^o hniq^o -D(q, w) used in the 
Einstein relation, Eq. ((43b|l . However, the diffusion re- 
sponse function D(q, w) was introduced in such a way 
that the limits lu —^ and g — > commute and we have 
only one definition for the static diffusion constant from 
the density response function 



D 



lim^^o V2x(q,w) 



q=0 



2 limq^o5t^x(q, ^)U=o 



2TmF 



lim,^og'$'^^(q,0) 



(47) 



The left equality in Eq. (|T7|) holds quite generally while 
the right one only at zero temperature. Note that repre- 
sentation ((46(1 holds only at zero temperature for small 
momenta and in the low- frequency limit. 

We can see from Eq. ((46(1 that the diffusion func- 
tion introduced in the quantum version of Pick's law 
reduces to the diffusion constant determining the long- 
range fluctuations of the electron-hole correlation func- 
tion and hence in the semiclassical limit also the diffu- 
sive behavior in the diffusion equation. However, only 
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the static diffusion constant D enters the denominator of 
the electron-hole correlation function. A generalization 
of the diffusion constant from the electron-hole correla- 
tion function to a frequency-dependent quantity via 

„^5H = ^V^$^>,c^)|^^„, (48) 

does not have a rigorous justification. It is nevertheless 
often used in the literature for a quantitative treatment of 
Anderson localization.^ Beware that the dynamical dif- 
fusion defined in this way does not fulfill the Einstein 
relation and hence deviates from the dynamical conduc- 
tivity. We show in the next section that in an exactly 
solvable mean-field limit the left-hand side of Eq. H48|l 
related to the conductivity via the Einstein relation dif- 
fers at finite frequencies from the conductivity obtained 
from the Kubo formula with the electrical current. 



VI. INFINITE-DIMENSIONAL MODEL: 
EXPLICIT EXACT SOLUTION 

We derived exact relations between density and cur- 
rent correlation functions for a lattice electron gas in a 



random potential. We now demonstrate the generally de- 
rived formulas explicitly on an exactly solvable limit of 
infinite spatial dimensions. This limit serves as a math- 
ematical tool for the definition of a mean-field theory 
not only for classical spin systems but also for itiner- 
ant disordered and interacting modelsi^^'^^'Si In case of 
the Anderson model of disordered electrons the mean- 
field theory, i. e., the limit of infinite spatial dimensions, 
equals the coherent potential approximation. We hence 
resort in this section to this solution. 

Equation determining the local self-energy in the 
coherent-potential approximation is Soven's equation 
that can be written as^^ 



G{z) = ([G-\z) + ^iz)~Vi\ (49a) 



where G{z) = N ^X^k^C'''^)- "^^"^ two-particle irre- 
ducible vertex then is 



A(z+,z_) 



G(z+)G(z_) 



1 + (I](z+) - Vi) G{z+) 1 + (S(z_) - V,) G(z_) 



I](z+)-S(z_) 
G(z+)-G(z_) 



AS 
AG 



(49b) 



The angular brackets were defined in Eq. The sec- 

ond equality in Eq. (|49bp is the Velicky-Ward identity. 

It is straightforward to find an explicit form of the 
two-particle Green function in the coherent-potential ap- 
proximation 



Gl,t^,(z+,z_;q) = G+(k)G-(k) [^(k-k') 

X{z+,z^)G+{k')G^(k') 



l-A(z+,z_)(G+G_) 



(50) 



where G±(k) = G(k ± q/2,z±) and (G+G_) = 
7V-iEkG+(k)G_(k). 

Using the two-particle Green function, Eq. 1)50(1 . we 
easily obtain explicit representations for the density- 
density and density-current correlation functions 



$(2;+,z_;q) 



$e(2:+,z_;q) 



(G+G_) 



l-A(z+,z_)(G+G_) 

(AqeG+G-) 
l-A(z+,z_)(G+G_) 



(51a) 
(51b) 



We now insert Eqs. H51|) together with the explicit rep- 
resentation for the two-particle Green function, Eq. (|50|l , 
to the right-hand side of Eq. H31a(l and find 



$,(z+,z_;q) - Az$(z+,z_;q) 



(Aqe - Az) 



A„G 



A„c-A2-AS 



l-A(z+,z_)(G+G_) 
(AG) - AS(G+G-) 
1 - A(z+,z_)(G+G_) ' 



(AG), (52) 



that is, continuity equation (|31a|l is fulfilled. Analogously 
we can explicitly verify Eqs. (|31b(l and H31c|l . Actu- 
ally, validity of continuity equations ((31|l follows from the 
Velicky-Ward identity, since the two-particle irreducible 
vertex X(z+,Z-) is momentum independent. In such a 
case the VoUhardt-Wolfle and the Velicky identities are 
equivalent. 

As next we show that the dc-conductivity in the 
coherent-potential approximation can be calculated using 
the diffusion constant D from Eq. To this purpose 
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we need to evaluate V^(G'+G_) at zero transfer frequency 
uj = 0. We have 

V^(G+G_) - iv,((G_VfeG+ - G+V,G_)> 

= i(G_v2G+ + G+V^G- - 2VfeG+VfeG_) 

= -(VfcG+VfcG_> (53) 

where in the last equality we used integration per parts 
in momentum space. We utilize this result and Eq. Hll|l 
to obtain the electron-hole correlation function 



$(q,Az) = — 



AG 



Az--- 



1 AS 



(54) 



2 AE - Az 



AG 



that in the low-energy limit reduces for real frequencies 
to 



<i>^«(q,c.) = 



27rnf 



(55) 



proving the Einstein relation for the dc-conductivity 
(jCPA p£ ^YiQ coherent-potential approximation. It is ev- 
ident from Eq. I|54(l that the frequency-dependent coef- 
ficient in the denominator of the electron-hole correla- 
tion function does not equal the dynamical diffusion from 
Eq. (|43b|l . To show this explicitly we use definition l|48|l 
and compare the result with the real-part of the conduc- 
tivity at zero temperature given by 



dE 
2Tr Lo 



- m^^iE, E + co; 0)] = a^^ico) + (c.) . (56a) 

We separated contributions to the conductivity from the 
electron-hole and electron-electron current-current cor- 
relation functions, It is useful to introduce 
also a simplified asymptotic form of the low-frequency 
conductivity 

SRa,„ (w) = -L [mi^ iEF^Lo,EF;0) 
Zn 

-m^^iEF-u;,EF;0)] (56b) 

where the integrand in the energy integral is replaced by 
its initial value. This representation is asymptotically 
exact in the limit uj —>■ for smooth current-current cor- 
relation functions ^aa near the Fermi energy Ef and 
carries the same frequency dependence as the diffusion 
constant D{uj) defined in Eq. 

For the numerical calculations we use a binary al- 
loy with two values of the random potential Vi — ±A 
weighted with probability x and 1 — x. Fig. ^ shows 
the two dynamical conductivities and the diffusion con- 
stant from Eq. (|48|l on a simple cubic lattice with pa- 
rameters X = 0.1, A — 0.6w for frequencies uj = 0.2w 
and UJ — OAw with w being the half-bandwidth. Note 




FIG. 1: Dynamical conductivities defined in three different 
manners, Eqs. H48|l . I|56|l for a binary alloy on a simple cubic 
lattice with the random potential Vt = ±0.6w, weighted with 
a; = 0.1 and 1 — x, respectively. Frequency is = 0.2'w in the 
upper pane and uj — OAw in the lower pane. Here w is the 
energy half-bandwidth. 



that the split-band value of the disorder strength is 
Ac « OAw. Although the three quantities differ for low 
frequencies only slightly inside the band, they behave 
differently near the band edges and for higher frequen- 
cies. In particular the diffusion constant shows anoma- 
lous behavior when the Fermi energy approaches a band 
edge. The compensating terms from the electron-electron 
(density-density) correlation function get relevant there. 
Anomalous behavior of the diffusion constant is more 
transparent for higher frequencies (lower pane). We can- 
not evidently rely on nFE>{uj) as a good approximation 
to the conductivity for finite frequencies except for Fermi 
energies deep inside the energy band. The smoothing im- 
pact of the integral over frequencies on the behavior of 
the dynamical conductivity gets clear from our numerical 
results. 

Recently a discussion was renewed about the propor- 
tion of contributions to the Kubo formula for the elec- 
trical conductivity from the electron-hole and electron- 
electron current-current correlation functions, and 
^BB. H56a|l . respectively.^^ In Fig. El these contri- 

butions are compared for the same setting of the binary 
alloy on a simple cubic lattice for w = 0. The contribu- 
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FIG. 2: Weight of contributions to the fuU conductivity from 
the electron-hole and electron-electron parts in the static case, 
to = 0. The setting is the same as in Fig. [U The lower pane 
shows the details of the satellite (split-off) band. 



FIG. 3: Electron-hole and electron-electron contributions to 
the conductivity at tj = OAw. The lower pane shows differ- 
ences between the conductivity obtained from the diffusion 
constant D and a^^. 



VII. DISCUSSION AND CONCLUSIONS 



tion from a dominates inside the band far from the 
band edges. Outside the central band and in the satel- 
lite impurity band the importance of the compensating 
effects of is evident. The situation worsens when we 
go over to the dynamical conductivity, Fig. 13 We can see 
that there is no region where the electron-hole contribu- 
tion would dominate or approximate the full conductivity 
reliably. It is important to note that conductivity a"^^ 
differs from the (dynamical) diffusion constant defined 
from the electron-hole density-density correlation func- 
tion, Eq. (|48|l . The latter contains even in the static limit 
both terms, and a^^. The difference is evident in 
Fig. 13 The conductivity calculated from D{oj) is much 
closer to the full electrical conductivity than a^^. With 
increasing frequency the role of the term increases. 
We hence cannot interchange contributions to the elec- 
trical conductivity from the electron-hole density-density 
correlation function <i>'^^(q, w) (its leading q-dependent 
term) and the conductivity o"^^(q, w). The former is 
generally better approximation than the latter except for 
Fermi energies near band edges and in band tails where 
both the approximations equally fail. 



In this paper we derived a number general relations 
between density-density, density-current and current- 
current correlation functions. With the aid of these re- 
lations we proved validity of an equation relating the 
electrical conductivity to the density response function, 
Eq. ^ , in the hydrodynamic limit of strongly disordered 
electron lattice systems and derived a dynamical general- 
ization of the Einstein relation between the conductivity 
and diffusion, Eq. . We found to what extent Eq. Q 
can be used in quantitative calculations and when it may 
be broken. We showed that the dynamical diffusion from 
the Einstein relation cannot be related to the diffusion 
pole beyond the static limit. We used only exact rea- 
soning and equalities derived directly from equations of 
motion for Green functions. We avoided extrapolations 
of semiclassical results or conclusions derived within equi- 
librium theory. 

The electrical conductivity can be related to the den- 
sity response function only if Ward identities are fulfilled. 
The same holds for the existence of a diffusion pole in 
the electron-hole correlation function with the static opti- 
cal conductivity governing its long-range spatial fluctua- 
tions. There are two, not fully equivalent, Ward identities 
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for noninteracting electrons in a random potential. Al- 
though they reflect conservation laws they can be proved 
only if perturbation expansion in the strength of the ran- 
dom potential for the irreducible two-particle vertex con- 
verges and results in an analytic function. 

In Sec. Hill we used Velicky-Ward identity Hill) to re- 
late mobility with the density response. This Ward iden- 
tity expresses conservation of probability in the state 
space and holds only for the homogeneous case, i. e., 
zero transfer momentum in the two-particle Green func- 
tion. To be of use in the calculation of the conductivity 
we had to assume analyticity of the hydrodynamic limit 
for all frequencies and to use the transfer momentum as 
an expansion parameter. VoUhardt-Wolfle-Ward identity 
(|25l) used in Sec. llVl to reveal the density response func- 
tion from the conductivity expresses charge conservation. 
It holds, unlike the Velicky-Ward identity, for arbitrary 
transfer momenta in the two-particle function, but it can 
be proved only via a perturbation (diagrammatic) expan- 
sion for the two-particle vertex function. Moreover, to 
warrant this identity we are not allowed to sum selected 
classes of relevant diagrams. A single relevant diagram 
for the two-particle irreducible vertex of order n produces 
n—1 irrelevant diagrams in the proof of the Ward identity. 
Selection of relevant diagrams may be, however, dictated 
by the existence of a singularity in the two-particle func- 
tion or by analyticity (causality) of the theory. There 
is hence no general guarantee that Ward identities and 
consequently the relation between the conductivity and 
the density response hold beyond the perturbative regime 
near the Anderson localization transition. 

When the conductivity can be related to the density 
response via Eq. Q we can derive a quantum version 
of the Einstein relation. To that purpose we introduced 
a diffusion response function relating the averaged cur- 
rent with the negative gradient of the charge density as 
a quantum generalization of Pick's law. Diffusion de- 
fined in this way enters the dynamical Einstein relation, 
Eq. H43|) , but it is not a priori clear whether this quantum 
diffusion reduces to the classical one in the semiclassi- 
cal limit. The solution of the classical diffusion equation 
pops up from a quantum theory via the diffusion pole and 
the long-range asymptotics of the electron-hole correla- 
tion function. We showed that the "classical" diffusion 
from the diffusion pole equals the quantum one from the 
Einstein relation and hence the conductivity only in the 
static limit. A dynamical diffusion constant defined from 
the electron-hole correlation function is no longer propor- 
tional to the conductivity beyond the static limit. This 
was demonstrated generally as well as explicitly on an 
exact solution in infinite spatial dimensions. 

To conclude, we showed that gauge invariance of an 



electron gas in a random potential responding to an elec- 
tromagnetic perturbation is guaranteed within linear- 
response theory if Ward identities Hll|) and (|25|l are 
fulfilled. In this case an (approximate) averaged two- 
particle Green function generates density and current re- 
sponse functions consistent in the hydrodynamic limit 
with relation Q). Or, what is more important, we can 
safely apply the approximate two-particle Green function 
only in either Kubo formula and use Eq. (0J to determine 
the other response function. We explicitly demonstrated 
that Ward identities are fulfilled only in the "diffusive" 
regime characterized by convergence of the perturbation 
expansion for the two-particle irreducible vertex and by 
analyticity of the hydrodynamic limit for all frequen- 
cies. We, however, know that in solutions with weakly 
or strongly localized electrons a pole in the two-particle 
irreducible vertex emerges and validity of the Ward iden- 
tities cannot be generally proved. One has to bear this 
aspect in mind when interpolating between the diffusive 
and localized regimes. 
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APPENDIX: PROOF AND VALIDITY OF THE 
VOLLHARDT-WOLFLE-WARD IDENTITY 

VoUhardt-Wolfle-Ward identity (|25|) is a primary tool 
for proving the continuity equations relating the cur- 
rent and density correlation functions. Understanding 
its proof and the domain of its validity is hence very im- 
portant. The original proof of VoUhardt and Wolflsl^. is 
restricted to real frequencies with different small imagi- 
nary parts only. However, it is valid for arbitrary com- 
plex frequencies. To demonstrate this we summarize the 
principal assumptions and steps of its proof from which 
we conclude on the domain of validity of Eq. H25() . 

First, we have to assume that the self-energy S is a 
functional of the full one-electron propagator G 

E(p,z) = S[G](p,z) (A.l) 

so that the perturbation expansion in the random poten- 
tial consists of irreducible diagrams only. We consider 
here only noninteracting electrons in a random poten- 
tial although the reasoning can be extended to corre- 
lated electrons as well. Perturbation expansion of the 
self-energy S(p, z) in the random potential reads 
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[(n+l)/2] „+l 



n=l k=l = i h,---,lk 



3k 



Ri,...,R„ + i V{ii,...,i^+i} 



— V e'^'i(P-Pi)e'^'2(pi-P2)e*R^„+i(P"-p)G'(p^^;z)...G(p„,z) (A.2) 

J n / J 



r 



where [x] is the integer part of a;, z„+i} denotes Only the momentum variables and momentum- 

a permutation of the indices {1, 2, . . . , n+ 1}, and the an- dependent functions are of importance. We denote I?„ 
gular brackets stand for cumulant averages defined from the sum of all diagrams with n internal fermionic lines 

with internal momenta pi, . . . ,p„. Then representation 
(jA.2|) can be simplified to 



(expl/),,=exp|f;l(T/"),| . (A. 



3) 



OO ^ 

5^(P, ^) = (^)a„ + E 1^ E E^'^"(^5P,Pl,---,Pn)G(pi,z)...G(p„,z) 



N' 

n—1 pi,...,p„ X>„ 



(A.4) 



We use notation G±(p) G(p ± q/2, z±) and S±(p) = S(p ± q/2, z±) and apply expansion (|A.4(I for the self-energy 
difference. We obtain 



OO 

AE(p) = s+(p)-E_(p) = ^— E^^"(^;P'Pi'---'P") 



n—1 pi,...,p„ Dn 

[G+(pi)...G+(p„)-G_(pi)...G_(p„)] . (A.5) 
The difference of the products of one-electron propagators can further be rewritten to a sum 



G+(pi) . . . G+(p„) - G_(pi) . . . G_(p„) = J2 G+(Pi) ■ • • G+(p,_i) 



X AG(p,)G_(p,+i)...G_(p„) (A.6) 



r 



where left (right) to the difference AG(pi) only G+(pj) 
(G_(pj)) appear. We sum all diagrams for the fixed dif- 
ference of the one-electron propagators, being now two- 
particle irreducible diagrams from the clcctron-holc chan- 
nel. We then come to a new representation 

OO n 

AE(p) - E E 1^ E ^ppT^' q) AG(p.) 



N 

n—1 i—1 pi 

(A.7) 

where App*-* (z+, z_; q) is a sum of two-particle irreducible 
diagrams with n internal one-electron lines of which i 
lines carry energy z+ and n — i lines energy z- . 

Last but very important step in the proof of identity 
H25|l is an assumption that removing the difference AG 



from Eq. (|A.7p and fixing the internal variable p; — p' 
does not change summability of the perturbation expan- 
sion. If so, we can write 

oc n 

App,(z+,Z-;q) = 5]5]A(;;^-')(z+,z_;q) . (A.8) 

n—1 i—1 

The equality holds if and only if the perturba- 
tion expansion in the random potential for the two- 
particle irreducible vertex App/ (z+, z_; q) converges 
point-wise for the chosen values of independent variables 
z+, z_, p, p', q. This is a rather severe restriction on ap- 
plicability of the Ward identity (|25|l . It says that the 
Ward identity is valid only in cases where no selective 
rules for sums of two-particle irreducible diagrams ap- 
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ply, i. e., all classes of diagrams are equally important. 
This happens when the irreducible two-particle vertex 
is regular (bounded). However, physically the most in- 
teresting situation occurs when due to backscatterings a 
Cooper pole appears in the irreducible vertex A. Thanks 
to this pole only a class of relevant (crossed electron- 
hole) diagrams determines the low-energy behavior of 
the two-particle vertex. In the asymptotic region of the 
Cooper pole in the electron-hole irreducible vertex, the 
Ward identity H25|) cannot be proved. In fact it gets vio- 
lated whenever we take into account a selected series of 
diagrams dictated by causality of the approximation or 
leading to a divergence and we have to consider a non- 
perturbative solution for the two-particle irreducible ver- 



tex. Each selected (relevant) diagram of order n from the 
formal expansion of the singular two-particle vertex gen- 
erates in the Ward identity n — 1 irrelevant diagrams to 
complete the sum from the right-hand side of Eq. ljA.6|l so 
that the self-energy difference can be represented via the 
vertex A. The irrelevant diagrams are neglected in ap- 
proximations with selected dominant classes of diagrams 
violating then identity (|25|l . Hence summability of the 
perturbation expansion for the two-particle irreducible 
vertex is an additive assumption in the proof of the Ward 
identity H25|l . There is no nonperturbative proof of this 
identity except for the homogeneous case, q = 0, where 
it is a consequence of the Velicky-Ward identity lfTT|l . 
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